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Some questions arising in the application of the thermal model to hadron production in heavy 
ion collisions are studied. We do so by applying the thermal model of hadron production to particle 
yields calculated by the microscopic transport model RQMD(version 2.3). We study the bias of 
incomplete information about the final hadronic state on the extraction of thermal parameters. It 
is found that the subset of particles measured typically in the experiments looks more thermal than 
the complete set of stable particles. The hadrons which show the largest deviations from thermal 
behavior in RQMD(version 2.3) are the multistrange baryons and antibaryons. We also looked at 
the influence of rapidity cuts on the extraction of thermal parameters and found that they lead to 
different thermal parameters and larger disagreement between the RQMD yields and the thermal 
model. 

PACS number(s): 25.75.Dw,24.10.Lx,24.10.Pa 



I. INTRODUCTION 

The hadron production in nuclear collisions at rela- 
tivistic and ultrarelativistic energies is under extensive 
experimental and theoretical investigation . One ques- 
tion of debate is whether hadron production is ther- 
mal, i.e., whether the hadron multiplicities for a given 
volume and energy density are distributed according to 
maximum entropy, or whether the particle multiplicity 
distribution reflects additional dynamical constraints. 

This question is not easy to answer. First the ther- 
mal model has various versions and depending on the 
additional input the answers on thermal parameters dif- 
fer |Q. Some of the differences in the various mod- 
els are the number of included resonance states the 
use of a canonical or grand canonical framework PB, 
the inclusion of hard core repulsion Q , the inclusion of 
strangeness suppression JjJ, the inclusion of other off- 
equilibrium parametrizations || , or the inclusion of finite 
size corrections in the density of states Q . 

The second important source of uncertainty in the 
judgment of thermal behavior is the incomplete measure- 
ment of the final hadronic state. Not all particles which 
are stable with respect to the strong interaction can be 
measured, and all measurements cover only a finite part 
of the complete momentum space. Some measurements 
cover nearly the full solid angle, e.g., NA49 Jl(J, in sym- 
metric collision systems, but some are restricted to only 
a small kinematic window |ll]-jl3| . 

The bias of the restricted information is the goal 
of our study here. We use the transport approach 
RQMD(version 2.3) |TJJ] for generating the full momen- 
tum space information on all stable particles. RQMD 



provides a microscopic description of heavy ion collisions 
which has been quite successful in predicting most of the 
observed features over a wide range of conditions and is 
therefore very well suited for our investigations. 



II. THERMAL MODEL 

The parameters of the thermal model may be under- 
stood as Lagrangian multipliers in extremizing the to- 
tal entropy under certain constraints. For a given point 
particle volume V we have as parameters the tempera- 
ture T (the Lagrangian multiplier to achieve a certain 
average energy density) and chemical potentials \iq (the 
Lagrangian multipliers to constrain the system on av- 
erage to a conserved total charge Q). In our case the 
conserved charges are baryon number B, strangeness S, 
and electrical charge Q. Instead of B, Q, and S one can 
equivalently conserve the net number of valence quarks 
for each flavor, (u — u), (d — d), (s — s), since on the rel- 
evant time scale only strong interactions are important. 
The relations between the chemical potentials in both 
descriptions are 

Mb = Mu + 2 ^d, 

US = Hd - Ms • (1) 

The corresponding fugacities are defined by 



A 



f = u,d,s 



(2) 



We know from pp collisions that strangeness is not in 
full chemical equilibrium with the nonstrange hadrons 
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|L5| . Therefore we allow for a deviation of strange par- 
ticles from complete chemical equilibrium. However, we 
force the strange particles to be in chemical equilibrium 
among each other pi. This requirement is equivalent to 
the introduction of the additional constraint that the to- 
tal number of strange valence quarks (s + s) be conserved 
on the average. The Lagrangian multiplier for this con- 
straint — expressed as a fugacity — is Jlq] 



7s = cxp(^ (s+5) /T) 



(3) 



7 S is a measure of the chemical saturation of strangeness 
with 7s < 1 meaning undersaturation, 7 S = 1 saturation, 
i.e., chemical equilibrium, and 7 S > 1 oversaturation. 

In the grand canonical description with the parameters 
introduced above the multiplicity Ni of a particle i with 
mass rrii and spin degeneracy Ji is given by 
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The (. . denote the corresponding valence quark con- 
tents of particle species i. The upper sign is for fermions, 
the lower sign for bosons. With the exception of pions we 
use for all particles the Boltzmann approximation; i.e., 
the sum over the index j is truncated at j = 1. 

For broad resonances (T > 1 MeV) we replace the 
S function in Eq. (^) by a Breit-Wigner distribution 
$Bw(rn 2 ) 
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where Ti is the width of the resonance i, mthrcs a thresh- 
old for the production of resonances, and 77i cut an upper 
cutoff. We use m t hrcs = ro» — 2r^ and m cu t = m, + 2Ti. 
The normalization a is determined by J dm 2 (*>bw( to2 ) 
1 . 

Equation (|^) gives the multiplicity of primary hadrons 
and resonances. The index for the particle species i 
runs over all hadrons of the Review of Particle Prop- 
erties with rtii < 1.7 GeV. The final hadron yields 
of stable hadrons include the decay contributions of the 
resonances. Unless otherwise stated we consider all these 
hadrons as stable which are stable with respect to strong 
decays. An exception is the <fi which is taken to be sta- 
ble in order to be in agreement with the special mode of 
RQMD used here (see below). The branching ratios are 
taken from [ fL7| , where missing or inaccurate entries have 
been completed by the requirement that all branching 
ratios for a resonance have to add to 1. 



The thermal model used here is in practical respects 
very similar to the one in ||. The main difference is 
the use of the grand canonical approach for all conserved 
charges in the study here. In addition there are small 
differences in the implementation of particle properties 
and decays where accurate numbers are missing in p7| . 
We use the following values for the assignment of hidden 
strangeness: 

(s + S) v = (s + s) r y = 1, 

(s + s)0 = (S + s)/ (980) = (* + s) /l( i420) = 2, 

(S + S) ^(1510) =(S + S) /^(1525) = ( s + S)^(1680) = 2 - (7) 

We show in Tables g and|| a comparison between the 
fits using the fully grand canonical approach and the 
partly canonical approach of ||. The two additional 
chemical potentials in our approach are determined by 
the same constraints as used in ||, i.e., strangeness neu- 
trality and electric charge conservation: 



SiN = 



-y 



BiN 



(8) 
(9) 



Z/A is the ratio of charge over mass number of the col- 
liding nuclei. We have the same number of degrees of 
freedom as in the canonical approach. We see from Ta- 
bles | and H that the two approaches are equivalent on 
the level of a few percent. 



III. THERMAL FITS TO 4tt MULTIPLICITIES 

We now turn to RQMD(version 2.3) as the source 
for our particle multiplicities. The RQMD transport 
model [[l4| is based on the excitation and fragmentation 
of colour strings followed by hadronic reinteractions. The 
version used [RQMD (version 2.3)] contains the color rope 
mechanism f25|| , an important source for the production 
of the strange hadrons which are of special interest in 
our study. It has been shown by comparison with ex- 
perimental data |26|,^7j that RQMD reproduces the main 
properties in hadron production. It is therefore a reliable 
model for generating the full momentum space informa- 
tion of all produced stable hadrons. In order to get also 
a <f> multiplicity we treated in this investigation the <p as 
a stable particle. 

We generated about 400 RQMD events each for S+S 
and Pb+Pb collisions with a centrality of 5% of crj°*j. 
The error on the input data is very important in the 
minimization of \ 2 . Hadron yields from RQMD have 
only a statistical error which can be made arbitrarily 
small if only enough events are generated. In order to 
put thermal fits to experimental data and thermal fits to 
RQMD yields on the same footing we give the multiplic- 
ities calculated by RQMD the same relative error as the 
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corresponding experimental value. The actual statistical 
errors from 400 RQMD events were in nearly all cases 
smaller than the corresponding experimental ones. We 
show in Table III a summary of relative errors used for 
RQMD multiplicities, cither taken directly from experi- 
ment or related to other entries of experimental relative 
errors. 

In Tables IV and |v| (rightmost column) the Air yields 
of the RQMD (version 2.3) simulation are shown as well 
as yields counted in a central rapidity window of varying 
size. The given multiplicities do not contain any decay 
contributions from weak or electromagnetic decays. The 
cited RQMD multiplicities differ in this respect from the 
experimental ones. In the thermal fits the corresponding 
changes have always been incorporated by the use of the 
appropriate decay channels. 

We performed different thermal fits to the RQMD Air 
multiplicities. First, we included in the thermal fit only 
the subset of particles as used in the thermal model de- 
scription of experimental data in H . The result is shown 
in Table VI. Compared to the thermal fit of experimen- 
tal data in |5j (see also Table |l]) we see that the thermal 
fit parameters of the RQMD analysis are close to the ex- 
perimental ones only in the case of S+S collisions. A 
significant discrepancy is seen for the Pb+Pb collisions. 
The extracted temperature for example differs by about 
30 MeV. We interpret this as an indication of larger dif- 
ferences in the hadron production between RQMD and 
the experiment. We do not try to investigate where the 
main differences between RQMD and the data lie, since 
this is, first of all, out of scope of this paper and, second, 
most of the experimental data are preliminary and thus 
it is too premature to speculate about physical reasons. 
The reader should further note that a direct compari- 
son of experimental yields and the here calculated yields 
from RQMD is not reasonable since the yields of K, A, 
and A contain no weak or electromagnetic contribution 
in the case of RQMD while such decays are partly in- 
cluded in the data. In addition, the experimental data of 
Pb+Pb collisions contain hadron yields and ratios which 
were measured in a limited acceptance and were extrap- 
olated to 4-7T. For all these reasons a thermal fit to the 
experimental Pb+Pb data would be premature. 

We would also like to draw the attention to the very 
low x 2 /DOF compared to the experimental fit. The 
small value of x 2 /DOF is partly a result of the used 
relative errors for RQMD. Since this error was taken 
from experiment, it contains the systematic error of the 
corresponding experiment. As a result of "perfect re- 
construction efficiency" in RQMD the error is overesti- 
mated. Using only the statistical error of RQMD would 
lead to a x 2 /DOF = 2.37/4 for S+S collisions and 
X 2 /DOF = 28.00/3 for Pb+Pb collisions. Neverthe- 
less, we keep the relative experimental error for RQMD 
yields, since we want to have the same relative weights 
of hadrons in a thermal fit as it is typical for an analysis 
of experimental data. We discuss the thermal character 
of RQMD yields in more detail below. 



We no w show in the rightmost column of Tables VII 



and VIII the resulting thermal parameters of a thermal fit 
including all stable particles from RQMD. The compar- 
ison between the RQMD yields and the thermal ones is 
given in Tables |^ and respectively. The most striking 
difference to the fit with the reduced number of hadrons is 
the increase in the x 2 /DOF by a factor of about 25, while 
the thermal parameters are basically the same. Only for 
7 S do we observe a small decrease with increasing number 
of stable particles. The result is understood by looking at 
the largest deviations of the thermal fit from the RQMD 
yields in Tables IV and [v|. In terms of standard devia- 
tions (s.d.) between the 4ir thermal fit and RQMD we 
have the largest deviations for (S+S, 29 s.d.; Pb+Pb, 
15 s.d.), Q (S+S, 9 s.d.; Pb+Pb, 8 s.d.), (S+S, 6.7 s.d.; 
Pb+Pb, 4.7 s.d.) and K+ (S+S 9 s.d.). Interestingly, 
also in the thermal fit of Ref. || to the Pb+Pb data 
the f2 and spoil the otherwise good fit, although their 
deviations go into the opposite direction. 

The above result shows that the thermal fit gets worse 
if all stable particles are considered, while a reduced num- 
ber of particle species suggest a better thermal behavior. 
The multistrange baryons are mostly responsible for the 
breakdown of the otherwise nice thermal behavior. This 
is illustrated in Figs, [j] and || where the result of the 
thermal fit with the reduced number of input particles 
of Table [vj is displayed. In addition to the fitted in- 
put particles (solid squares) the resulting other multi- 
plicities of these thermal fits are shown as open circles. 
For both collision systems the fi, Cl, and S show the 
largest deviations. The x 2 /DOF of the thermal model 
with the parameters of Table [vl applied to all stable 
multiplicities of Tables IV and jyL respectively, leads to 
X 2 /DOF = 224/21 for S+S_and x 2 /DOF = 191/21 for 
Pb+Pb, while excluding £1, f2, 2°, S - , S°, and S~ leads 
to x 2 /DOF = 8.5/15 for S+S and to x 2 /DOF = 6.9/15 
for Pb+Pb. 

We see here the breakdown of the concept of relative 
strangeness equilibrium [Q which is violated by various 
strange particle species. The thermal fit has the ten- 
dency to underestimate the meson yields, and the ther- 
mal \S\ = 1 and \S\ = 2 (anti)baryons are roughly in 
agreement with RQMD, while the 15*1 = 3 (anti)baryons 
are overpredicted by the thermal model with optimized 
parameters relative to RQMD. The decrease of 7„ from 
the thermal fit with a reduced number of hadrons com- 
pared to the fit with all stable hadrons is a result of the 
above-mentioned undersaturation of the RQMD multi- 
strange hadrons. This can nicely be seen in Figs. ^ and 

An explanation for the above findings was offered in 
[g7| where the unexpectedly small inverse slope Q of 
the pr spectrum was investigated within the RQMD 
model. It was found that the f2 decouples very early from 
the rest of the fireball. It is therefore very unlikely that 
the Q comes into chemical equilibrium with the rest of 
the hadrons. The same may also partly be true for the 
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A word of caution is necessary here. We see in Figs, [l] 
and H that the Q and f2 have the largest deviations from 
the thermal fit. However, one must note that the physics 
of the n(Q) production in RQMD is not treated on the 
same high quality level as the other particles. Since the 
O was of minor physical relevance before — this has 
changed with the high precision measurement of f2 ra- 
tios |3Cf | and slopes |28| — a precise implementation of 
the production physics in RQMD was neglected. The 
used strangeness suppression factors and decuplet-octet 
suppression factors result already in e + e~ collisions in 

and f2 yields which are a factor of 2 too low. So far 
no f2(51) production via resonance states has been incor- 
porated. The above-mentioned shortcomings reduce the 
f2(f2) yields. Therefore an improvement of the fi physics 
within RQMD will lead to smaller deviations from the 
thermal fit than seen now. 

We conclude from this section that a relatively small 
subset of hadrons looks more thermal than the complete 
set of stable particles. Thermal behavior should there- 
fore always be tested taking as many hadron species into 
account as possible. Nevertheless, we saw only minor 
changes in the thermal parameters between the thermal 
fit with a reduced number of hadrons and the thermal 
fit with all stable hadrons. Therefore, the experimen- 
tally accessible subset of hadrons is suited to determine 
the thermal parameters temperature, chemical potentials 
and the fireball point particle volume. ~f s , however, is 
sensitive (at least in the case of RQMD) to the multi- 
strange (anti)baryons and should therefore be extracted 
by considering the full set of particles. 



IV. THERMAL FITS INCLUDING RAPIDITY 
CUTS 



Particle yields or ratios are very often only accessible 
in a limited kinematic range. While an extrapolation in 
Pt can be performed rather reliably due to the observed 
nearly exponential behavior of the mx spectra, an ex- 
trapolation in rapidity is much more difficult because dif- 
ferent hadron species are known to have quite differently 
shaped rapidity spectra. A thermal fit in case of kine- 
matic cuts needs model assumptions about the momen- 
tum distributions of the particles. The simplest assump- 
tion one can make is the assumption of Bjorken scaling 
in rapidity |pl|| . However, as a result of the much too low 
beam energy, Bjorken scaling has not been seen in RQMD 
or in the experimental data at SPS energies. Taking, 
however, the analyzed rapidity window to be very nar- 
row around midrapidity the Bjorken scaling assumption 
might be approximately valid. In ]32| ] it is shown that 
within the hydrodynamic model the extraction of ther- 
mal parameters from yields in a limited rapidity range 
works quite reasonably. We like to investigate the sys- 
tematic error of rapidity cuts in a similar way with the 
help of RQMD events. 



For this investigation we first integrate the particle mo- 
mentum distribution over various rapidity intervals (see, 
e.g., Table 0). Then we assume Bjorken scaling around 
midrapidity in the sense that these yields or ratios are fit- 
ted by Eq. (^) without correction due to the finite rapid- 
ity interval. Therefore the fitting procedure is the same 
as for the full momentum space yields, but the volume 
has to be considered as the sum of comoving subvolumes. 



We show in Tables IV and VII t he re sults for central 
S+S collisions and in Tables |Vj and VIII the correspond- 
ing results for Pb+Pb collision. In S+S collisions we 
see a clear decrease of the temperature with decreasing 
width of the rapidity window which is accompanied by 
an increase in 7 S and a decrease in X u ^d- The decrease 
of X u .d is a result of the baryon hole around midrapid- 
ity in S+S collisions p3fl . The increase in 7 S reflects the 
fact that RQMD, like the experiment, produces strange 
quarks more efficiently in the central region. On a weaker 
scale the same tendencies are seen for Pb+Pb collisions. 

An interpretation of this result is done in connection 
with a similar study of S+S collisions within the hydro- 
dynamic model p2[ . In this study much smaller devi- 
ations between a global thermal fit and particle yields 
resulting from integration over different rapidity win- 
dows were found. In order to understand the difference 
one has to realize that in RQMD the physics, especially 
the physics of strange particle production, in the cen- 
tral region is different from the fragmentation regions. 
Therefore the change here in thermal parameters with 
the analyzed rapidity window resembles mostly a change 
in physics, while in the case of the hydrodynamical model 
by construction the pure kinematic bias was investigated. 
While a purely kinematic influence increases the temper- 
ature in the central region [ j32| the differences in physics 
decrease the central temperature compared to a Air anal- 
ysis. It was further found |52| that for systems with large 
longitudinal flow the artificial change of thermal parame- 
ters is small. Therefore the change of thermal parameters 
with the rapidity window in the case of RQMD is clearly 
dominated by the physics component. 

Thus we interpret the large change of thermal param- 
eters with different rapidity window size as a result of 
no global equilibrium in these collisions and the devia- 
tions from a global equilibrium is much larger in S+S 
collisions compared to Pb+Pb collisions. However, there 
might still be local thermal and chemical equilibrium in 
the experiments as well as in the RQMD model. For 
a study of local equilibrium in RQMD we refer to J|J] 
and in UrQMD to |35| . We only want here to emphasize 
that rapidity cuts can lead to considerable changes in the 
thermal parameters adjusted to a collision system. The 
change is the larger the further away the system is from 
global equilibrium. In addition we see that the x 2 /DOF 
is best for 47r integrated input data; as the rapidity in- 
terval is narrowed, we see a fluctuating but increasing 

x 2 /dof. 
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V. CONCLUSIONS 



We have used the RQMD model to address some ques- 
tions arising in the application of the thermal model to 
particle production in ultrarelativistic heavy ion colli- 
sions. We have performed a thermal fit to RQMD mul- 
tiplicities, where the errors in the RQMD yields have 
been adjusted to be of the same order as typical for ex- 
periments. First we found that for an input consisting 
of only a small number of hadron species — typical for 
existing experiments — the RQMD yields can be nicely 
reproduced by a thermal model. This suggests at a first 
glance global thermal behavior. However, comparing the 
thermal model with the full set of stable particles leads to 
significant deviations. The largest discrepancy from ther- 
mal behavior is found for the multistrange hadrons. The 
usual small sample of particle species still gives nearly 
the same thermal parameters as the fit with the full set 
of stable hadrons. Experimentally the 47r data are very 
scare and extrapolation is always required when a ther- 
mal fit is performed. Our results imply that one must 
be careful when using the thermal model to analyze the 
experimental data since the constrain is global, not local 
{§■ 

It is difficult to assess whether the complete set of 47T 
multiplicities including the multistrange (anti) particles 
from RQMD shows still "good" thermal behavior. A 
judgment on the basis of the absolute value of x 2 /DOF 
is not reasonable because x 2 /DOF depends strongly on 
the assumed error of the multiplicities and, given certain 
existing deviations of the RQMD output from thermal 
behavior, becomes arbitrarily large as the statistical er- 
ror on the yields becomes smaller and smaller by gener- 
ating more and more events. A more useful criterion for 
judging the quality of the thermal fits may be given by 
the average deviation of the best thermal fit (with the 
relative errors from experiment) from the RQMD mean 
multiplicity. It is of the order of 20% for the fit in Tables 
[V and |v|, of order 8% for a fit excluding the £1 and f2, 
and of order 4% for a fit excluding fi, all S states and 
the corresponding antiparticles. Given that the absolute 
yields of the different hadrons species vary over three or- 
ders of magnitude these numbers suggest that RQMD 
yields are roughly in agreement with a thermal model 
interpretation. However, it must be shown — similar 
to the studies for the kinetic freeze-out in RQMD 63] 
or as done for the UrQMD model in j35fl — that the ex- 
tracted chemical freeze-out parameters are meaningful in 
the sense that the microscopic state of RQMD at chem- 
ical freeze-out is in agreement with chemical and kinetic 
equilibrium and that the resulting microscopic thermal 
parameters are the same as for the analysis presented 
here. 

In a very recent publication the result of a thermal fit 
to UrQMD midrapidity ratios was presented |37j . In the 
case of Pb + Pb 158A GeV collisions an agreement was 
found between the UrQMD multiplicities and the ther- 



mal model on the same quality level as presented here 
but with slightly different parameters for temperature 
140 MeV H vs 154 MeV here (see Table [VHl| left col- 
umn) and fi B /T = 1.50 [0 vs 1.57 here. The detailed 
microscopic study in |37j revealed that two fundamen- 
tal assumptions of the "thermal model," i.e., universal 
freeze-out and global equilibrium, are violated; from this 
the authors J57| concluded that thermal fits to particle 
ratios are meaningless. We, on the other hand, would like 
to reserve judgement on the meaningfulncss of thermal 
fits to particle multiplicities until we understand why a 
thermal fit to particle yields works much better than the 
large violation of common freeze-out (27jj3lj] and global 
equilibrium |57 35 would suggest. 

We also investigated the influence of rapidity cuts on 
the extracted thermal parameters. There is a change of 
thermal parameters with decreasing size of the analyzed 
rapidity window. This change is dominantly due to non- 
global equilibrium effects in these collisions. Thus one 
has to keep in mind that there is a nonnegligible bias on 
the extracted thermal parameters due to cuts in rapidity. 
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TABLE I. Comparison between fitted and measured hadron abundances and ratios. Both thermal model calculations include 



the weak/electromagnetic decay of S°, 3°, 3 , Q. and the corresponding antiparticles. 



Hadron 


Measured 


Ref. 


Fitted in 


5| Grand canonical fit 


S+S collisions 


h- o 


98±3 




92.63 


94.21 


K+ 


12.5±0.4 




12.68 


12.65 


K~ 


6.9±0.4 




7.611 


7.27 


K° 


10.5±1.7 




9.939 


9.72 


A( b ) 
A( b ) 


9.4±1.0 




7.692 


8.30 


2.2±0.4 




1.474 


1.633 


p-p n 
p( d ) 


21.2±1.3 




21.49 


22.47 


1.15±0.4 


■21 


2.092 


2.044 


Pb+Pb collisions 


Net baryon 


372±10 


mmm 


375.7 


374.3 


h- (») 


680±50 




650.2 


661.0 


K° 


68±10 


IIlI 


58.27 


58.22 


4> 


5.4±0.7 


■ill 


5.759 


5.660 


p-p ( c ) 


155±20 


lul 


155.3 


152.8 


K+/K- 
A/A( b ) 


1.8±0.1 


mmm 


1.652 


1.699 


0.2±0.04 


l£l 


0.188 


0.195 



a - Defined as ir~ + K~ + p. 
b - Includes feeding from 3. 

c - Measured with the '+' - '-' method, in this case limited 
rapidity acceptance (0.2-5.8) to exclude spectators. 
d - Measured in a restricted rapidity interval and extrapolated, 
by assuming that p has the same rapidity distribution as the A. 



TABLE II. Hadron gas model fitted parameters. The fit itself is shown in Table M. The entries denoted by an asterisk are 



derived quantities while the others are fit parameters used by MINUIT. 



parameter 




S+S 






Pb+Pb 


ref. [5 




grand canonical 


ref. [5] 


grand canonical 


T (MeV) 


182.4±9.2 




187.7±1.4 


192.6±8.1 


188.8+1.1 


y T 3 e -0.7GcV/T 


3.51±0.15 




3.42* 


24.3±1.6 


24.4* 


V (fm 3 ) 


206* 




166±11 


992* 


1134+76 


Is 


0.732±0.038 




0.731±0.033 


0.616±0.043 


0.620+0.047 




1.248±0.074 




1.251* 


1.207±0.071 


1.258* 


An 






1.516± 0.021 




1.469+ 0.018 


A d 






1.518± 0.021 




1.547+ 0.022 


\s 






1.069± 0.005 




1.071+ 0.005 


X7D0F 


17.1/4 




12.0/4 


3.99/3 


2.33/3 
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TABLE III. List of relative errors of 4-zr multiplicities and ratios used for RQMD multiplicities. 







s+s 






Pb+Pb 


particle 


rel. error in % 




source 


rel. error in % 


source 



p-p 

net baryon 
K+/K" 
A/A 



3.1 
6.1 



p 


D.l 


p 


O4.o 


11 


D.l 


11 


O^l.O 


■7T + 
71 


3 1 


7T 


'■i 1 


_o 


3 1 


K+ 

IV 


'4 9 


IV 


. O 


K° 


1 9 


rh 




A 


in 

1U.D 


A 


18.2 


£+ 


10.6 


£+ 


18.2 


E° 


10.6 


£° 


18.2 


E~ 


10.6 


£" 


18.2 




10.0 




15.0 




10.0 




15.0 




23.0 


Q- 


34.5 



as p-p 



as Pb+Pb 



7.4 
12.9 
2.7 
5.6 
20.0 
12.9 
34.8 
12.9 
34.8 
7.4 
7.4 
7.4 
3.2 
5.8 
14.7 
13.0 
10.6 
18.2 
10.6 
18.2 
10.6 
18.2 
10.6 
18.2 
10.0 
15.0 
10.0 
15.0 
23.0 
34.5 



( a ) We assume that the relative error is the same as the error on multiplicity 
per event in the kinematic window of WA97 for the most central events |2S 

( b ) The relative error of multiplicities of multistrange antibaryons is typically 
50% larger than the relative error of the corresponding baryons |29|. 
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TABLE IV. RQMD(version 2.3) multiplicities of 200 A GeV/c S+S collisions. Trigger cross section a trig < 5%a gcom (b < 0.85 
fm). 





lyl <o.5 




\y\ <i o 




\v\ <i-5 




\y\ <2.o 




47T 




RQMD 


therm. 


RQMD 


therm. 


RQMD 


therm. 


RQMD 


therm. 


RQMD 


therm. 


V 


3.13±0.19 


3.297 


6.95±0.42 


7.304 


11.65±0.71 


11.954 


16.47±1.00 


16.624 


24.24±1.47 


23.552 


V 


0.50±0.17 


0.342 


0.99±0.34 


0.975 


1.25±0.43 


1.122 


1.36±0.47 


1.187 


1.40±0.48 


1.301 


n 


3.20±0.19 


3.292 


7.23±0.44 


7.296 


11.79±0.71 


11.943 


16.60±1.01 


16.611 


24.23±1.47 


23.536 


n 


0.49±0.17 


0.340 


0.89±0.30 


0.972 


1.12±0.39 


1.119 


1.23±0.42 


1.183 


1.27±0.44 


1.297 


7T + 


18.42±0.57 


18.719 


35.78±1.10 


36.959 


50.36±1.56 


51.527 


60.89±1.88 


62.484 


72.54±2.24 


74.863 


7T~ 


18.80±0.58 


18.732 


36.10±1.11 


36.983 


50.56±1.56 


51.564 


60.18±1.86 


62.532 


72.69±2.25 


74.919 


7T° 


19.18±0.59 


20.562 


37.15±1.15 


40.267 


52.31±1.62 


56.141 


63.38±1.96 


67.837 


75.31±2.33 


80.891 


K + 


3.34±0.10 


3.136 


6.40±0.20 


5.477 


8.92±0.28 


8.018 


10.51±0.33 


9.384 


11.45±0.36 


10.354 


K~ 


2.10±0.12 


2.152 


3.86±0.22 


3.702 


5.16±0.29 


4.988 


5.91±0.34 


5.469 


6.26±0.36 


5.567 


K° s 

4> 


2.70±0.43 


2.623 


5.08±0.82 


4.563 


6.98±1.13 


6.467 


8.11±1.31 


7.386 


8.76±1.41 


7.921 


0.42±0.054 


0.309 


0.81±0.10 


0.519 


1.06±0.13 


0.754 


1.17±0.15 


0.811 


1.19±0.15 


0.790 


A 


1.05±0.11 


0.838 


2.21±0.23 


1.663 


3.32±0.35 


2.688 


4.18±0.44 


3.416 


4.72±0.50 


4.200 


A 


n o K -1-C\ n A c 
U.Z0±U.U40 


0. 131 


n A A _l_rt non 
U.44±U.U8U 


0.343 


n c c _l_n i A 
U.o6±0.10 


0.429 


n cn±n i n 
0.6U±U. 1U 


0.447 


n a i _l_rt 11 
U.D1 ±U. 1 1 


0.469 


v+ 


U.OU±U.Ooo 


0.329 


u.yo±u. ru 


0.650 


i /in _l_n i c 
1.4y±U. 10 


1.050 


1.84±U.19 


1.333 


011 _l_rt oo 


1 .636 


Zj ' 


U. 1U±U.018 


0.051 


U.17±0.U3U 


0.134 


0.2o±U.041 


0. 168 


U.25±0.U45 


0.175 


0.25±U.045 


0. 183 


v° 


U.4o±U.U4o 


0.325 


n 014.(1 fifif; 


0.643 


l.zyztu. lo 


1 .039 


i .oy±u. 10 


1.320 


1.1 1 ZtU. 10 


1.620 


V D 


U.11±U.UZU 


0.050 


n i n_l_n nQ /i 


0.133 


r» n A -1-C\ rt A o 
U.z4±U.U4o 


0. 166 


n oc_l_n n .1 *7 
U.26±0.U4 ( 


0.173 


0.2o±U.04 / 


0.181 


E~ 


0.47±0.049 


0.319 


0.94±0.099 


0.633 


1.39±0.14 


1.023 


1.74±0.18 


1.300 


1.92±0.20 


1.597 


E~ 


0.11±0.020 


0.049 


0.19±0.034 


0.130 


0.24±0.043 


0.163 


0.26±0.047 


0.170 


0.26±0.047 


0.179 


— 


0.11±0.011 


0.135 


0.25±0.025 


0.237 


0.36±0.036 


0.378 


0.43±0.043 


0.439 


0.44±0.044 


0.466 




0.07±0.010 


0.032 


0.12±0.018 


0.075 


0.15±0.022 


0.102 


0.16±0.024 


0.105 


0.16±0.024 


0.105 




0.14±0.014 


0.133 


0.26±0.026 


0.233 


0.36±0.036 


0.372 


0.43±0.043 


0.432 


0.44±0.044 


0.459 




0.06±0.009 


0.031 


0.11±0.016 


0.074 


0.13±0.019 


0.101 


0.15±0.022 


0.103 


0.15±0.022 


0.103 


si- 


0.01±0.002 


0.016 


0.01±0.002 


0.028 


0.02±0.004 


0.044 


0.02±0.004 


0.047 


0.02±0.004 


0.044 


a- 


0.01±0.003 


0.005 


0.01±0.003 


0.014 


0.01±0.003 


0.020 


0.01±0.003 


0.021 


0.01±0.003 


0.020 



TABLE V. RQMD(version 2.3) multiplicities of 158A GeV/c Pb+Pb collisions. Trigger cross section ctrig < 5% a gC0 m (b < 3 
fm). 





\y\ <o. 




\y\ <i 


.0 


\y\ <i- 






.0 


4tt 




RQMD 


therm. 


RQMD 


therm. 


RQMD 


therm. 


RQMD 


therm. 


RQMD 


therm. 


V 


27.06±3.49 


32.321 


55.08±7.10 


64.884 


82.51±10.6 


96.019 


105.00±13.5 


120.752 


131.40±17.0 


146.992 


V 


1.68±0.58 


2.295 


3. Mil. 09 


4.011 


3.92±1.36 


5.184 


4.37±1.52 


5.746 


4.58±1.59 


6.348 


n 


28.89±3.72 


33.513 


58.00±7.48 


67.476 


85.82±11.1 


100.184 


109.35±14.1 


126.259 


138.83±17.9 


154.054 


n 


1.77±0.61 


2.203 


3.25±1.13 


3.839 


4.05±1.40 


4.946 


4.45±1.54 


5.471 


4.66±1.62 


6.030 


7T + 


144.16±10.7 


156.050 


271.73±20.1 


291.882 


369.94±27.4 


397.705 


435.64±32.2 


470.565 


500.05±37.0 


534.792 


7T~ 


151.11±11.2 


163.344 


284.80±21.1 


306.556 


389.28±28.8 


419.270 


460.63±34.1 


497.486 


531.42±39.3 


567.035 


n° 


152.56±11.3 


176.860 


287.54±21.3 


331.296 


391.99±29.0 


450.767 


461.97±34.2 


532.048 


531.37±39.3 


601.890 


K+ 


31.91±1.02 


30.007 


60.16±1.92 


56.971 


81.12±2.59 


76.911 


93.69±2.99 


88.895 


100.33±3.21 


95.780 


K~ 


18.76±1.08 


19.155 


34.72±2.01 


35.180 


45.55±2.64 


45.836 


51.51±2.98 


51.683 


54.28±3.14 


54.358 


K°s 


25.59±3.76 


24.488 


47.92±7.04 


45.934 


63.73±9.36 


61.239 


72.75±10.7 


70.173 


77.30±11.4 


75.0 


<t> 


4.55±0.59 


3.179 


8.19±1.06 


5.999 


10.74±1.39 


7.853 


11.82±1.53 


8.654 


12.08±1.57 


8.681 


A 


10.79±1.14 


8.936 


20.63±2.18 


17.769 


29.34±3.11 


25.326 


35.38±3.75 


30.367 


38.09±4.03 


34.206 


A 


1.24±0.22 


1.019 


2.06±0.37 


1.829 


2.62±0.47 


2.367 


2.90±0.52 


2.569 


2.99±0.54 


2.696 


S+ 


4.30±0.45 


3.412 


8.49±0.89 


6.770 


12.20±1.29 


9.623 


14.48±1.53 


11.518 


15.32±1.62 


12.946 


S+ 


0.55±0.10 


0.411 


0.88±0.16 


0.740 


1.12±0.20 


0.960 


1.25±0.22 


1.043 


1.27±0.23 


1.096 


s° 


4.21±0.44 


3.465 


8.50±0.90 


6.890 


12.07±1.27 


9.817 


14.35±1.52 


11.769 


15.38±1.63 


13.250 


s° 


0.47±0.085 


0.395 


0.83±0.15 


0.709 


1.07±0.19 


0.918 


1.17±0.21 


0.996 


1.20±0.21 


1.045 


E~ 


4.45±0.47 


3.496 


8.81±0.93 


6.967 


12.35±1.30 


9.952 


14.77±1.56 


11.949 


15.76±1.67 


13.478 


E~ 


0.49±0.089 


0.377 


0.86±0.15 


0.676 


1.14±0.20 


0.873 


1.27±0.23 


0.945 


1.28±0.23 


0.991 


— 


1.57±0.15 


1.532 


3.01±0.30 


3.008 


4.13±0.41 


4.114 


4.66±0.46 


4.696 


4.78±0.47 


4.880 




0.42±0.063 


0.295 


0.72±0.10 


0.545 


0.95±0.14 


0.707 


1.04±0.15 


0.752 


1.05±0.15 


0.751 




1.53±0.15 


1.528 


3.04±0.30 


3.006 


4.16±0.41 


4.118 


4.71±0.47 


4.706 


4.84±0.48 


4.896 




0.41±0.061 


0.284 


0.69±0.10 


0.524 


0.92±0.13 


0.679 


1.00±0.15 


0.722 


1.00±0.15 


0.720 


Sl~ 


0.10±0.023 


0.203 


0.20±0.046 


0.396 


0.26±0.059 


0.525 


0.28±0.066 


0.571 


0.28±0.066 


0.552 


n~ 


0.03±0.010 


0.064 


0.06±0.020 


0.122 


0.08±0.027 


0.160 


0.08±0.027 


0.167 


0.08±0.027 


0.159 
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TABLE VI. Thermal parameters of fitting RQMD yields and ratios in S+S (left) and Pb+Pb (right) collisions. The same 
particles and ratios were taken like in the fit of ||. The input data set was extracted from the entries in Tables [Fv]and 
respectively. The relative error is as in Table |. 



S+S Pb+Pb 



particle 


RQMD 


thermal 


particle 


RQMD 


thermal 


hr 


82.98+2.57 


83.77 


h~ 


612.5+45.32 


598.87 


K° 


8.76+1.42 


8.74 


J *■ s 


77.31+11.36 


81.71 


p-p 


22.84+1.39 


22.22 


p-p 


126.82+16.36 


126.50 


K+ 


11.45+0.37 


11.53 




12.0+1.57 


11.88 


K~ 


6.26+0.36 


6.02 


K+/K- 


1.847+0.104 


1.832 


A 


4.72+0.50 


4.82 


A/A 


0.0785+0.016 


0.080 


A 


0.61+0.11 


0.597 


net baryon 


346.6+9.36 


347.19 


P 


1.4+0.49 


1.38 








X7DOF 




0.83/4 






0.29/3 


T(MeV) 




179.4+7.5 






156.7+8.6 


V(fm 3 ) 




195+65 






3008+1780 


Is 




0.769+0.045 






1.129+0.17 


A (( 




1.609+0.026 






1.677+0.044 


A d 




1.610+0.026 






1.761+0.065 


A s 




1.102+0.032 






1.208+0.059 



TABLE VII. Thermal parameters of fitting RQMD multiplicities in various rapidity intervals for S+S collisions (trigger 
cross section crtrig < 5%a gcom ). The fit itself is shown in Table |lV| 





\y\ <o.5 


Ivl <i-o 


Ivl <1.5 


\y\ <2.o 


4-k 


T(MeV) 


153.3+1.3 


166.6+3.9 


167.5+2.7 


168.5+2.7 


172.9+3.0 


V(fm 3 ) 


147+8 


173+29 


227+26 


261+30 


260+32 


7s 


0.905+0.027 


0.768+0.022 


0.793+0.019 


0.748+0.018 


0.670+0.016 




1.460+0.017 


1.400+0.016 


1.485+0.015 


1.555+0.017 


1.623+0.019 


Ad 


1.461+0.017 


1.401+0.016 


1.487+0.015 


1.556+0.017 


1.625+0.019 


A s 


1.188+0.007 


1.124+0.010 


1.136+0.010 


1.145+0.011 


1.138+0.013 


XVDOF 


115.32/21 


152.20/21 


98.51/21 


111.91/21 


94.72/21 



TABLE VIII. Thermal parameters of fitting RQMD multiplicities in various rapidity intervals for Pb+Pb collisions (trigger 
cross section crtrig < 5%a gcom ). The fit itself is shown in Table 





| V | <0.5 


Ivl <i-o 


Ivl <1.5 


\y\ <2.o 


4% 


T(MeV) 


152.8+4.8 


153.3+3.3 


154.4+2.9 


154.4+2.7 


155.4+2.6 


V(fm 3 ) 


1223+277 


2215+347 


2868+390 


3386+427 


3706+460 


7s 


1.019+0.046 


1.024+0.037 


1.000+0.033 


0.966+0.031 


0.900+0.029 


A ( ( 


1.538+0.032 


1.573+0.024 


1.608+0.022 


1.641+0.022 


1.667+0.022 


A d 


1.595+0.041 


1.635+0.031 


1.676+0.028 


1.714+0.028 


1.745+0.028 


A s 


1.211+0.017 


1.216+0.014 


1.219+0.014 


1.227+0.014 


1.230+0.014 


XVDOF 


79.62/21 


66.17/21 


67.38/21 


67.61/21 


55.61/21 
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FIG. 1. Thermal fit of RQMD 4-k yields of S+S collisions. The solid squares denote the input yields of the reduced fit of 
Table VI. The open circles show the resulting other 4-7T multiplicities of stable particles not included in the fit. 
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FIG. 2. Thermal fit of RQMD An yields of Pb+Pb collisions. The solid squares denote the input yields of the reduced fit of 
Table VI. The open circles show the resulting other 4tt multiplicities of stable particles not included in the fit. 
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